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The optimal power performance of a first principle quan- 
tum heat engine model shows friction-like phenomena when 
the internal fluid Hamiltonian does not commute with the 
external control field. The model is based on interacting two- 
level-systems where the external magnetic field serves as a 
control variable. 

I. INTRODUCTION 

It is well established that the performance of work- 
ing heat engines are limited by intrinsic unavoidable ir- 
reversibilities. Maximum power is obtained at the ex- 
pense of efficiency where the reversible point of maxi- 
mum efficiency has zero power. This principle has been 
clearly illustrated by the endoreversible model of Curzon 
and Ahlborn [Q and summarized by Salamon et. al. [^. 
Two additional unavoidable sources of loss are heat leaks 
which practically eliminate the maximum efficiency adia- 
batic operation, and internal friction which restricts fast 
operating cycles. 

Is this universal performance limitation of heat engines 
macroscopic or microscopic? Though the common image 
of heat engines is of large macroscopic devices, micro- 
scopic models based on first principle quantum mechanics 
are limited by the Carnot efficiency [Q, and show a re- 
markable resemblance to their macroscopic analogs when 
the engines produce finite power |^]. 

In previous studies, both discrete |^-^ and continuous 
quantum models PJTo|-[T2t have been scrutinized, which 
are analogous, respectively, to four-stroke-engines and 
turbines. The present study examines a discrete four- 
stroke quantum engine, comparing it to an engine subject 
to phenomenological internal friction. It will be demon- 
strated that the quantum engines inability to control si- 
multaneously the external and internal portions of the 
working fluid Hamiltonian is its source of friction. 

II. BASIC CONSTRUCTION 

All working heat engines and all refrigerators operate 
on the same principle. The engine manipulates the en- 
ergy flow between three reservoirs, a hot, cold and power 
reservoir, either to extract power from the temperature 
difference or to pump heat from the cold to the hot reser- 
voir at the expense of external power. 

The present model of a discrete quantum heat engine 
is composed of a cycle of operation constructed from two 
adiabats and two isotherms similar to a Otto cycle. The 



quantum dynamics are generated by external fields in the 
adiabats and by heat fiows from hot and cold reservoirs 
in the isochores. The working medium is modeled as a 
gas of interacting particles with the Hamiltonian: 

H = Hext + Hint • (1) 

Hext = ^ Hi is the sum of single particle Hamilto- 
nians, lu — Lu{t) is the time dependent external control 
field, and Hi„t represents the inter-particle interaction. 

The change in time of an operator A during the adia- 
batic/isochoreal branches is described as: 

A = z[H,A] + ^ + J^niA) , (2) 

Cd represents the Liouville dissipative generator on the 
isochore in contact with either the hot or cold bath (units 
of h = 1). Replacing A by H in Eq. (^) leads to the 
power invested in or extracted from the adiabats:, 

V = (Hi) , (3) 

i 

where (Hi) is the expectation value of the single particle 
Hamiltonian. The heat flow is extracted from the energy 
balance on the isochores |^,^ : 

Qh/c - (A/c(h)) . (4) 

For simplicity, the single particle Hamiltonian is chosen 
as a two-level-system (TLS): lli = a^. The interaction 
term is restricted to coupling of pairs of spin atoms. As a 
result, the state of the working medium is described by an 
ensemble of pairs of two-level-systems represented by the 
density operator p, defined in the tensor product space 
of the individual two TLS systems. Expectation values 
are obtained by the usual definition (A) = tr{Ap}. The 
external Hamiltonian then becomes: 

Hext = u(^al(g,P +P (g,al^ (5) 

and the external field is chosen to be in the z direction. 
The interaction Hamiltonian is chosen as: 

Hint = J {al®al-al®al) . (6) 

J scales the strength of the interaction. When J 
the model approaches the previously studied friction- 
less model |7j. The inter-particle interaction term, Eq. 

, defines a correlation energy between the two sin- 
gle particle spins in the x and y direction. As a result, 
[Hext, Hint] 7^ 0, since the external Hamiltonian is po- 
larized in the z direction. 



1 



III. DYNAMICS OF THE WORKING MEDIUM 

The dynamics generated by the Hamihonian is com- 
pletely determined by the algebra of commutation rela- 
tions of the set of 16 operators spanning the Hilbert space 
of the combined system. Due to symmetry, the commu- 
tation algebra decomposes into subsets of operators with 
closed commutation relations. The set is generated by 
commutation relations between the operators composing 
the Hamiltonian: 



Bi 
B2 



- 1 



(g) + (g) 

- 2 - 1 „ - 2 



(7) 



On the isochores, the system is in contact with a ther- 
mal bath which eventually will lead the working fluid to 
thermal equilibrium with temperature T: 



-/3H 



Peq = 



Z 



(11) 



with /3 = l/fc&T and Z = tr{e^'^^}. The dynamics 
generated by the system-bath interaction is described by 
the dissipative Liouville operator CD^ which in Lindblad 
form becomes [[l3|: 



E 



FiXFl - -(FiF^X + XFiF^ 



The commutation relation: [Bi, B2] — 4iB3 leads to the 
definition of B3 which closes the set i.e. the set of oper- 
ators Bi, 62,63 forms a closed subset of the Lie algebra 
of the combined system. The Hamiltonian expressed in 
terms of the set of operators becomes: H = ojBi + JB2 

The commutation relations of the set of Bj operators 
are isomorphic to the angular momentum commutation 
relations when the transformation -jBi Jj is applied. 
This similarity can be exploited to express the expecta- 
tion values in a Cartesian three dimensional space, where 
the external field is in the Bi or z direction, the correla- 
tions in the B2 or x direction and B3 is in the y direction. 

The closed set of operators Bi is sufficient to follow the 
changes in energy and to obtain the power consumption. 
Using Eq. (||) and the commutation relations of the set 
of Bj operators, the Heisenberg equation of motion for 
this set becomes: 



d 
di 



These equations can be written in matrix form for the 
expectations bi = tr{'Bi,p}: 
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Since the matrix A{t) is time dependent, the propagation 
is broken into N short time segments At where NAt — t, 
and is solved numerically. The matrix A is diagonal- 
ized for each time step, assuming that toit) is constant 
within the time period At. The corresponding eigenval- 
ues become:— 4in, and Ai^l, where ft = \jlP^+u? . 
The short time propagator for the adiabats from time t 
tot + At: 



Ua{t,At) = e 



_ pA(t)At 



w,/(l-c) 

— Js 

n 




(10) 



where c = cos(4f2At) and s ^ sm{4ilAt). 



where Fj are operators from the Hilbert space of the sys- 
tem. The operators Fj which control the approach to 
thermal equilibrium become the transition operators be- 
tween the energy eigenstates. 

Substituting the Bi operators into Eq. (O) one gets: 



rBi + 2^(fci-fct) 
rB2 + |i(fci-fct) 
r B3 



(13) 



where F — fc 4 + fc |, and the coefficients k | and k t 
obey detailed balance !^■ = e~^'^^, with the bath 



temperature (3 = -ff-j-- The set of Bi operators and the 

identity operator I form a closed set to the application 
of the dissipative operator £d. 

The relaxation to equilibrium is accompanied by loss 
of phase. Additional dephasing can be caused by elastic 
bath fluctuations which modulate the systems frequen- 
cies. This pure dephasing conserves the systems energy 
£d''(H) — 0. It is obtained by inserting the Hamiltonian 
in Lindblad's form ( Eq. ( |l2| ) ) as one of the operators 
Fi: 



Ld\^) = 7[H,[H,X]: 



(14) 



Equations of motion for the set of Bi operators on the 
isochore are obtained from Eqs. (|l|) (|l|) and (|): 



where: 



dt 



r + 167 J2 
—1670-17 

-4J 



Bh 



— 16"fJu! 

-r -I- 167^2 

4w 



(15) 
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(16) 



The solution of Eq. ( p^ for the isochores becomes: 

h{t) = UriHO) -hS'i) + h^i (17) 
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where b^'^ 



jtC and 



-(r-i67a^)t 
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where X — exp{—16jfl^t), 
sin(4m). 



Js 

4 I ' (18) 



= cos{Ailt) and s 



IV. THE CYCLE OF OPERATION 

Fig. 1^ illustrates the cycle of operation on the plane of 
the variables to, the external control, and the projection 
of the polarization on the energy axis E/il. A different 
view is displayed in Fig. |^ showing the cycle trajectory 
on the volume defined by the set of "polarization" coor- 
dinates (Bi). 




FIG. 1. The heat engine with J 7^ in the u , E/Q plane. 
Th is the hot bath temperature, th is the time allocation when 
in contact with the hot bath. Tc and Tc represent the temper- 
ature and time allocation for the cold bath. Tf,a represents the 
time allocation for compression (field change from LJb to uja) 
and Tab for expansion. The rectangular box which includes 
point E on the hot isochore where the system is in equilib- 
rium with Th, and point F (equilibrium with Tc ) on the cold 
isochore, represents the cycle of maximum work. This cycle 
spends an infinite amount of time on all branches. The work 
is the area of the rectangle Wmax = (fit — f2a)-(-§^ — ^). The 
optimal power cycle is emphasized by the green shading. The 
time allocations for this cycle are th = I.705,r(,a — 0.4369 
,Tc = 1. 76597321 and Tat = 0.4953. The cycle depicted in 
purple is characterized by very short time allocations on the 
adiabats (r^a ~ Tab = 0.00025). The common engine parame- 
ters are: uJa = f .794, ujb = 4.239 Th = 2.5, Tc = 0.5, J = 0.6. 
Fc — Th = l with units where h — 1 and kb — 1. 

The cycle starts at point A where the system is in con- 
tact with the hot bath at temperature Th- The system 
accumulates heat for a period until it reaches point B. 



Point E is the equilibrium point of the scaled energy E/^ 
at the bath temperature Th with external field strength 
cjf). Adiabatic compression from to Ua follows the tra- 
jectory from point B to point C for a time duration r^a 
with a constant co. The system is in contact with the 
cold bath from point C to point D for a time duration 
Tc- The cold bath equilibrium point, at temperature Tc 
at LUa, is F. The cycle becomes closed by a compression 
stage from point D back to point A. 




FIG. 2. The optimal power heat engine cycle corresponding 
to Fig. I. The Bi coordinate represents the individual parti- 
cle polarization, and the B2 and B3 coordinates represent the 
interparticle correlation. The inner cycle with blue isochores 
is subject to strong dephasing (7^ = 0.15 and 7c = 0.5). 

For long time duration on the adiabats the cycle of op- 
eration is restricted to the Bi, B2 plane. For fast motion 
on the adiabats the system cannot follow the instanta- 
neous change in the direction of the Hamiltonian rotat- 
ing on the Bi, B2 plane. As a result the expectation of 
the B3 operator increases starting a precession type mo- 
tion around the temporary direction of the Hamiltonian 
H = u)(t)'Bi + JB2. This can be seen clearly in Fig. || in 
the trajectory from point D to point A. The precession 
motion continues on the isochore where the Hamiltonian 
becomes constant (point A to point B). In addition, due 
to dephasing, the amplitude of the precession motion is 
damped. Part of the dephasing is caused by the energy 
equilibration with the bath. When pure dephasing is 
added the precession motion is damped almost instantly 
(Cf. fig 1^). The motion out of the Bi, B2 plane causes 
a bending upward in the adiabats as seen in Fig. 0. This 
bending causes additional work on the adiabats which is 
then dissipated on the isochores. This causes a reduction 
in efficiency from rjmax = 1 ^ ^a/^b ~ 0.5581 which is 
reached at infinite cycle times to a lower value at maxi- 
mum power r]Pmax = 0.495 (rjcamot = 0.8). The cycle of 
the engine is completely determined by the external con- 
trol parameters uja,uJb,Th,Tc, and the time allocations: 
Tab- Independent of the initial condition, the 
engine settles to a limit cycle after a few revolutions with 
the preset sequence of isochores and adiabats. 
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The optimization objective is the power of the engine 
which is the total work per cycle W divided by the cy- 
cle period r. Work is obtained only on the adiabats 
and is calculated by integrating the instantaneous power 
Eq. (HI) for the adiabat duration: Wba = Jq*"^ Vdt = 

The optimization analysis starts by setting the external 
parameters as the extreme field values iVb, (^a and the hot 
and cold bath temperatures at Th, Tc- The performance 
of the engine therefore will be determined by the time 
allocated to the different segments. By setting the total 
cycle time t = Tta + Tc + Tab + tu the optimization is 
carried out by partitioning the time between the adiabats 
and isochores. This splits the allocated time between the 
hot and cold bath isochores, and splits the allocated time 
between the compression and expansion adiabats. 

In the limiting case of no internal coupling J = 0, the 
current model is identical to the noted frictionless one 
0,1). In this frictionless case the optimal power time al- 
location on the isochores becomes: ThTh = FcTc and zero 
time allocation on the adiabats. For J 7^ the time allo- 
cations changes considerably. Two limiting cases emerge, 
the slow limit where most of the cycle time is allocated 
to the adiabats, and the fast or sudden limit where most 
of the time is allocated to the isochores. 

Due to the precession motion, the cycle operation is 
noisy. A small change in a parameter can considerably 
alter the limit cycle and thus the performance of the en- 
gine. Additional pure dephasing damps the noise of the 
engine (Cf. Fig. |^) while the overall performance is only 
shghtly altered. 

The global power optimum was sought by both a con- 
jugate gradient method and by a random search scruti- 
nizing local maxima. The optimal power as a function 
of the total cycle time is shown in Fig. ^ for different J 
values. It is clear that the power has a clear maximum 
with respect to the cycle period t. The optimal value de- 
creases and cycle time increases with increasing J. The 
maximum power output as a function of J is shown in 
the insert together with the analogous friction case. 

Despite the large local fluctuations with respect to time 
allocations the optimal power performance of the quan- 
tum engine shows a remarkable overall similarity to the 
performance of an engine subject to phenomenological 
friction as studied in Ref. Q . One expects friction to op- 
pose the fast motion on the adiabats, therefore the extra 
power invested has to be independent of the sign of the 
change in the control field. This means that to the lowest 
order it has to be proportional to the square of the time 
derivative i.e. Vfric — o^'jP'. The accumulated extra 
work is then dissipated on the cold isochore. The engine 
subject to friction shows performance curves and opti- 
mal time allocations which are very close to the present 
first principle quantum model. For this case the origin of 
lost power is the inability of the "polarization" vector to 



follow adiabatically the instantaneous Hamiltonian. The 
resulting precession motion on the isochores then leads to 
additional dissipation on the cold isochore. To the lowest 
order, the additional power should scale as which ex- 
plains the observed linear relation of J with the friction 
parameter a. 
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FIG. 3. Power as a function of cycle time r for optimal 
time allocation on the four branches for three values of the 
parameter J. The underlying dotted lines are the optimized 
power output of an engine with phenomenological internal 
friction (in black) . A linear relation between J and the friction 
parameter a fits the data. The insert shows the Maximum 
power as a function of J (in red) and a x 30. 

To conclude, we have found that a quantum heat en- 
gine with a working fluid which is not completely con- 
trollable by the external field shows performance charac- 
teristics which can be mapped into a heat engine subject 
to phenomenological friction. 
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